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Abstract 

An excitonic method proper to study conjugated oligomers and 
polymers is described and its applicability tested on the ground and first 
excited states of trans-polyacetylene, taken as a model. From the Pariser-Parr- 
Pople Hamiltonian, we derive an effective Hamiltonian based on a local 
description of the polymer in term of monomers; the relevant electronic 
configurations are build on a small number of pertinent local excitations. 
The intuitive and simple microscopic physical picture given by our model 
supplement recent results, such as the Rice and Garstein ones. Depending of 
the parameters, the linear absorption appears dominated by an intense 
excitonic peak. 
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I- Introduction 



71 electrons confer to conjugated polymers attractive electronic 
properties which are associated to the traditional elastic ones of saturated 
polymers. For instance, the optical absorption thresholds of these 
compounds are in the visible or near UV range. Some of them possess a 
strong and fast non-linear response. Last, some compounds present an 
insulator-metal transition under doping. During the past twenty years it 
became clear that the electron-phonon and the electron-electron 
interactions are both important to correctly describe the lowest electronic 
excitations of n electrons [1]. Moreover the long-range part of the electron- 
electron interaction term has somehow to be considered, in order to 
properly describe, for instance, the excitonic states clearly observed in 
polydiacetylenes (PDA). The simplest and more studied model Hamiltonian 
which includes these specific interaction terms is the so-called Pariser-Parr- 
Pople Hamiltonian [2] 



H PPP ~ 52 t n,n+^ a na a n+\,a + a n+\o a n o ) + 
n,a 

— V„ m ( c£„a„„ )(a + a 

2 n,m* no no „ ma , 



(1) 



ma 2 ' 



n,m,a,a 



In this expression a na (a n(J ) are the creation (annihilation) operators of an 
electron in site n with spin o. The first term describes the kinetic energy of 
the 7C electrons and their interactions with the atomic cores, expressed using 
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first neighbours hopping integrals. The second term is the Coulomb 
repulsion between % electrons localized on the sites n and m. Various 
parametrizations of the coulomb term exist in the literature [1]. Here in 
view to compare with the results of Yu et ah, we adopt their parametrization 
[3] 

V 

V nm = Uiin = m and V nm = :if n * m. (2) 

\n — m\ 

The electron-phonon interaction is treated semi-classically by introducing a 
linear dependence over the bond lengths in the hopping integrals. In the 
simplest case, the polyacetylene (PA), an usually adopted parametrization 
was introduced by Su, Schrieffer and Heeger (SSH) [4, 5] 

tn,n+1 =to+a{Un-Un+l) (3) 

where a is the electron-phonon coupling constant and u n is the 
displacement co-ordinate of the nucleus n along the molecular axis. The 
bond alternation observed in PA is easily reproduced by setting u n = (-1) n u, 
yielding a pair of integrals f3 d and f3 s , the hopping integrals associated with 
the double (1.35 A) and the single bonds (1.45 A) respectively. 

Hppp is invariant under electron-hole transformation. Consequently, 
the eigenstates of this Hamiltonian can be classified following the 
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electron-hole symmetry classes i.e. the so called (+) and (-) classes. This 
symmetry allows to simplify the calculations. 

Overview of the methods employed for the study of H PPP . With the 
parameter values relevant for conjugated polymers, the coulomb terms are 
approximately of the same order of magnitude than the kinetic ones. In this 
case, the study of the PPP Hamiltonian becomes a very difficult problem. 
Only the shortest polyenes (with the number of double bonds N, N < 8) can 
be properly studied by the usual quantum chemistry methods [2]. For the 
thermodynamic limit additional, and often drastic, approximations must be 
applied. The mean field theories [4, 5] gave a simple physical picture of 
conjugated polymers but, since electronic correlations are discarded, the 
results are always quantitatively incorrect and even, in some cases, 
qualitatively wrong [1]. A direct refinement of these methods are the 
excitonic ones, introduced for instance by Abe et al. [6]; using the valence 
and conduction bands obtained by mean field theory, they perform a 
configuration interaction including every monoexcited Slater's 
determinant. For a given parameters range, exciton states may split out of 
the electron-hole pairs continuum. With such a method the optical 
absorption spectrum of PDA, in which exciton states have been 
experimentally observed, is rather well reproduced [6]. However, electronic 
correlations are still missing; a resulting failure can be found in the loss of 
size consistency when biexcited states are introduced [7]. Recently Yu et al. [3] 
combined the former theory [6] with a band calculation using the projection 
technique of Becker and Fulde [8, 9] to include correlation effects. They study 
the first excitations of the polyparaphenylenevinylene (PPV) considered as 
an effective linear chain. Good agreement is found with experimentally 
known energies of the first singlet and triplet excitons as well as with the 
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threshold of the conduction band. Moreover, inclusion of biexcited states 
brings no size consistency problems contrary to the previous methods. 
However, with the procedure of renormalization of Yu et ah, the two 
particles of hole-electron pair are independently renormalized. This 
procedure becomes questionable when the two quasi particles are 
constrained to stay close to each other as in the exciton states of conjugated 
polymers. We will examine more precisely this point later. 

Molecular crystal approach. On the other hand, a more intuitive way is 
to consider the polymer chain as an one-dimensional molecular crystal of 
monomers. This approach was first proposed by Simpson [10], refined by 
many authors [11-13] and recently used by Rice and Garstein [14-16] on 
polyparaphenylene (PPP) and PPV. Grounds can be found in the prevalent 
simple bond character of the intermonomer bonding. 71 electrons are then 
preferentially confined on the monomers according to their specific 
topology [17]. In these simple models photoexcitations of the polymers are 
derived from the local excitations of the monomers and from charge 
transfer excitations between monomers. Analytical results can then be 
obtained and a good agreement is found with the absorption spectrum of 
PPP, after adjustment of the parameter values [16]. The same authors 
considered also the PPV as a PPP with breaking of the electron-hole 
symmetry; again, a good agreement is found for the linear absorption [18]. 
However, the Rice and Gartstein's model (RG) is less microscopic than the 
one described by H PPP . Moreover the field of application of the RGs model 
is limited to the study of optical absorption spectrum although 
Mukhopadhyay and co-workers have recently employed a similar 
molecular exciton method in view to describe the spin-charge cross-over in 
dimerized chain [19]. 
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Scope and plan. In this paper, we present a new method of calculation 
for the ground state and the first excited states, lending a particular attention 
to the linear absorption. This new method aims to bridge the gap between 
the traditional quantum chemistry or solid state physics methods based on 
Hppp [2, 3] and the simpler molecular exciton methods [14]. Purposely, the 
electronic configurations of the polymer are built on the basis of localized 
self-consistent orbitals of the monomers. This specific choice of basis is a 
natural one for molecular exciton methods. Then from H PPP we build an 
effective Hamiltonian by selecting the particular subspaces of electronic 
configurations relevant to the ground state and the primary excitations. 
With this procedure, the excited states are eventually obtained with the 
same formal expression as in the RG model. However, with our method, 
the empirical parameters of RG model are expressed in function of the ones 
of Hppp and the physical understanding of the states is ameliorated as we 
will see below. Also, our calculations save size-consistency as the ones of ref 
[3] but remain much more simple. For simplicity we develop here this new 
method for the simplest conjugated polymer, the trans-PA in the neutral 
state; extensions to more complex compounds as PPP and PPV are 
straightforward and will be presented elsewhere. 

The outline of this paper is the following. In section II we introduce 
the self-consistent orbitals of the monomer (ethylene) and the interaction 
terms corresponding to this special choice of one-electron functions. We 
present then the configuration subspace which will be used as model space 
for the ground state (section III). An approximate diagonalization into this 
subspace will be given in section IV. Finally we describe in section V the first 
excited states of (+) electron-hole symmetry - relevant for the linear 
absorption - and we will discuss our results. 
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II- Description of the polymer from the self-consistent orbitals of the 
monomers 

In first approximation conjugated polymers are quasi-one 
dimensional compounds with a carbon backbone characterised by several 
types of bonds. For instance, an alternation of the lengths of the different 
bonds, double and single, occurs in PA, the simplest conjugated polymer. 

o o 

The double bond (1.35A) is shorter than the single one (1.45A) and, due to 
the linear dependence of the hopping terms, |/^|>|/? s . |. Because of this 

general feature, and in the spirit of the molecular exciton methods, we 
choose the self-consistent orbitals of the monomers as one-electron base 
functions. This choice of representation for the polymer is the first step 
towards an approximate description of the ground state and of the first 
excited states. The critical parameter which controls the relevance of this 
particular representation is the bond alternation parameter z = ^ where a is 

the projection of the average C-C spacing onto the chain axis and £ is the 

electronic coherence length defined as t = ^° a. When z = (£ — > °°), the 

2au 

chain would not be dimerized, with equal lengths for every bond. This case 
is the more inappropriate to be dealt with our model. Beneficially, larger the 
order alternation parameter is, and more our particular basis will constitute 
a good starting point for the study of the % electronic properties. When z = 1 
= a) the chain is totally dimerized, no electron transfer occurs from a 
monomer to another and our special choice is the best one. 
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For the ethylene simple symmetry considerations give the expressions 
of the self consistent orbitals. We associate with these two orbitals the 
following creation (destruction) operators for the ethylene n 




na — 



42 



1 



( a2na + a(2n+i )a) ) associated 



with the bonding orbital ; 




a(2n+i)a)) associated 



with the antibonding orbital. 

A polyene with N double bonds is then represented by N two levels 
systems (Fig. 1). 

On this basis, the energy terms of H PPP , can be reorganised in three 
different classes (Fig. 2): 

(i) The kinetic term which allows an electron or a hole to hop from a 
monomer to one of its first neighbours. This term is given by (3 S multiplied 
by a constant which depends on the topology of the system (and of the 
orbital symmetries). This constant is equal to 1/2 for the poly acetylene. 

U-V 

(ii) The intramonomer coulomb term which introduces 

electronic correlation inside the monomer, by coupling fundamental 
configuration of a monomer with the doubly excited one. 
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(iii) The intermonomer coulomb dipolar terms are of two types: 

- First the transition dipole-transition dipole interaction terms 
between two monomers distant by r monomer units. One effect of these 
terms is to allow the interaction of intramonomer monoexcitations. The 
corresponding expression is 



- Second the transition dipole-permanent dipole interaction which is 



These dipolar interaction terms decrease very rapidly on r and can be 
considered as important only for r=l. 

In their model Hamiltonian, Rice and Garstein have introduced some 
semi-empirical parameters which are not considered in the traditional 
mean field theory [14-16] : 



r(') = -^( V l,2r ~ 2V l2r+1 + V 1>2r+2 ) 



(4) 



expressed by 




(5) 
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- the "correlation energy gap" U RG which is the energy required in the 
limit j3 s — > to dissociate a local monomer excitation into a well separated 
electron and hole; 

- the dipole-dipole interaction V RG between local excitations on 
neighbouring monomers. 

With our formulation of the Pariser-Parr-Pople Hamiltonian, for the 
case of the polyacetylene these semi-empirical parameters acquire a more 
microscopic meaning. The identification of V RG with r(1) is immediate but 
the interpretation of the U RG term is less direct and will be given in 
section V. 

Ill- Generative local electronic configurations - Collective excitations 

Generative local configurations. In our localized picture of the 
polyacetylene, electronic configurations will appear as combinations of local 
electronic configurations. To generate a tractable model, a subset of the 
Hilbert space has to be selected. A simple way to do it is to retain only some 
relevant local electronic configurations. In this chapter, we will illustrate 
this selection procedure in building the ground state subspace. For excited 
states it will be done in section V. 

The ground state traditionally adopted by the molecular excitonic 
method is [12, 19] 
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N 



1 ) = J"| 5+ t | Vacuum ) (6) 



where jVacwwra) is the vacuum state, i.e. the state without any n electron. 
The energy of this state will serve as reference in the following. The state \0) 
only includes one type of Local Configuration (LC). We will name it F-LC 
and the corresponding creation operator is 



K = b+ aKi (7) 



At this level of approximation, the monomers are considered as 
independent and each monomer possesses 2 n electrons which are described 
in the mean field approximation. With this very simple picture the 
dynamics of the n electrons would be poorly described. In particular, the 
conjugation phenomenon proper to the 71 systems is not reproduced and the 
electronic correlations are not introduced inside each monomer. We then 
improve this description by introducing two electronic local configurations 
able to interact directly with the state \0) by one of the characteristic 
interaction terms described above. 



1- The D-LC in which the monomer n is doubly excited ; this LC is 
associated with the creation operator : 
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D, 



— A + -t>A + i B a.B i 

n\ n\ n \ ni 



(8) 



This LC interacts with the state \0) by 



U-V 
2 



and introduces intramonomer 



electronic correlation; the corresponding energy is E D = 4j3 d . 

2- The Tcf-LC in which one electron is transferred from monomer n 
(or n + 1) to monomer n + 1 (or n); the transfers on the left and on the right 
are combined with a (-) sign. This particular linear combination belongs to 
the class of electron-hole symmetry noted (-) which is the ground state 
symmetry class. This LC is associated with the creation operator: 



This LC directly interacts with the state \0) through f$ s and introduces short 
range intermonomer charge fluctuations. This effect is intended to correct 
the too strong localization of the electrons on the double bonds associated 
with \0); the energy of this LC is E(1) = 2j3 d + V + A(1) where A(r) is the 
attractive interaction between the hole and the electron at a distance r 




(9) 



Mr) = - 1 -{v h2r + 2V l!2r+1 + V 1i2r+2 



) (rX>) 



(10) 
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These two kinds of local excited configurations bring about a local 
improvement of the 71 electrons dynamics. Therefore in order to improve 
the dynamics of the whole tc electronic system, it is necessary to consider 
electronic configurations containing a number of 7c 7 ~-LCs and D-LCs, able 
to introduce simultaneously intermonomer electron derealization and 
intramonomer coulomb correlation in any location of the chain. Indeed, we 
have to build electronic configurations by combining the three LCs F-LC, D- 
LC and Tq-LC, which, in this sense, will be called generative local 
configurations (GLC). 

A given electronic configuration involves n t TcJ-LQ located on sites 
labelled {x(i)} (i = 1,..,n t ) and n D D-LC located on sites labelled {y(j)} 
(j = 1,..,n D ); the remaining monomers are represented by F-LC 

\x(1),...,x(nj;y(1 ),..., y(n D )) = Tc+ (1) ...Tc+ (nt) D; (1) ...D; (nD) \0). (11) 

More precisely, a Tcj -LC configuration extends on two next neighbour sites; 
the label x(i) is defined for the left site. In (11) spatial overlap between GLCs 
are forbidden, so that x(i) # x( j ),x( j )±1 V(z,y) and y( j ) * x(i),x(i) + 1 

V(i,y)- 

The energy of these electronic configurations is given by the following 
expression: 
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E({x(i);y(j)}) = n t E(1) + n D E D + W({x(i)}). 



(12) 



Because of the long range part of the coulomb potential, the energy of an 
electronic configuration depends on the relative positions of the n t Tq-LC 
through the term ^({jcf i)}). This dependence is in fact very smooth so that 

in practice we will neglect it. We can write : 



E{{x( i);y( j)}) « E(n t ,n D ) = n t E( 1) + n D E D (13) 



Collective excitations. From the electronic configurations (11), we build 
collective excitations which are characterised by {n t ,n D }, respectively the 

numbers of Tcf-LCs and D-LCs. n t and n D are not really independent each 
other; indeed, once the n t Tcf-LCs have been localized on the chain 
(n t = 0,...,N / 2), it is obviously not possible to place a D-LC on the sites 
labelled {x(i)} or {x(i) + 1}. Consequently n D = 0,...,(N-2n t ). For the 
thermodynamic limit the relevant collective excitations are : 



\{n D ,n t })= . 1 S\\x(1),...,x(n t );y(1),...,y(n D )). (14) 

lx (i£i(n,n 

\y(1),...,y(n D )l 



Where N( n D ,n t ) is the number of electronic configurations with n t Tc-f -LCs 
located on sites {j f }and n D D-LCs on sites {x D }. In these expressions the 
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summation is over the whole electronic configurations that it is possible to 
perform. This summation is symbolised by { }. 

In the expression (14) the determination of the normalisation constant 
is purely a problem of enumeration. First let us deal n t TcJ -LQs on the N 
double bonds of the polymer. Every Tq-LC occupy two neighbour 
monomer sites, on the both sides of a single bond. The problem is 
equivalent to place n t Tcf-LCs on (N-l) single bonds but with the additional 
constraint to avoid their overlap. So, it is necessary to introduce between 
two 7cy~-LCs a forbidden zone which extents at least over one single bond. 
With this topological constraint n t independent LCs have to be placed on an 
effective chain of N effective sites. N is equal to the number of single bonds 
(N-l) minus the number of forbidden bonds (n t -1); so N = N-n t and the 
configuration number for n t 7cy~-LCs is : 

^(n t Tc1-LC) = cl = C N n - n '. (15) 

Note that this enumeration can be also obtained by recursion. 

It is simpler to introduce the on-site D-LCs. The problem is to place 
independently n D D-LCs (n D = ,N '- 2n t ) on the remaining (N-2n f ) 
monomers. This is a classical result 

*an D D;(N-2n t )) = C»- 2n >. (16) 
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Following (15) and (16) we finally complete the description of the collective 
excitation 

N(n D ,n t ) = C»-*C»; 2 *. (17) 
IV- Approximate ground state 

Generally speaking, determining the exact ground state of H PPP is a 

very difficult task which can only be performed numerically for finite 

compounds with less than six double bonds. In the polymer limit, various 

more or less drastic approximations have been proposed such as the 

simplified ground state of the molecular exciton methods (6) [12, 13, 19], the 

traditional mean field ground state [6] or the sophisticated Giitzwiller 

variational solution of the Peierls-Hubbard Hamiltonian [20]. In this section 

an approximate ground state is built by diagonalizing H PPP in the Hilbert's 
subspace spanned by the collective excitations |{«o,« f }). We believe that this 

approximation grasps the essential features of the ground state in view to 
describe the linear spectroscopic properties of conjugated polymers. 

The collective excitations expressed by (14) and (17) interact each other 
through two distinct interaction terms : 
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l tW> n D) = ({ n D> n t + 1 }\ H PPp\{ n D> n t}} 



(18-a) 



7 Z)(«D '> n t) = {{ n D + 1 > n t} \ H PPP | { n D > n t }) ( 18 " b ) 

the intermonomer delocalization interaction term and the intramonomer 
correlation interaction term respectively. The dipolar terms which have 
been described above are small enough to be neglected in this approximate 
treatment of the ground state. 

Intramonomer electronic correlation. First we consider the 

configuration subspaces which are spanned by the collective excitations 
|{ft£>,« f }), such as n D varies from to (N-2n t ); each subspace is 

characterised by a specific value of n t . Only the second interaction term acts 

inside a given subspace : 

I D {n D ;n t ) = ^{n D + l){N-2n t -n D ) 1 ^- (19) 

where n D = ,...,( N '- 2n t ) , that allows an easy diagonalization of H PPP in 
these subspaces. Indeed, from a particular electronic configuration with n t 
localized Tcf-LCs \x(1 ),..., x(n t )), we can independently introduce 
intramonomer electronic correlation on each of the (N-2n t ) remaining 
double bonds [21]. We then obtain for each monomer two states |-) and |+) 
associated with the energies £_ and £ + respectively 
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e ± =2p d ±yi6pl + (U-V) 2 



(20) 



If the state \0 n ) in which the particular monomer n is in its own ground 
state is introduced, the states |-) and |+) are written as 



-) = a \0 n ) + bD+ \0 n 

(U-V) l — p 

where a = ' ' and b = i1 - a 2 . (21) 

\ + ) = aD:\0\-b\0\ J4e* + (U-Vf 



From the electronic configuration \x(1 ),..., x(n t )), we obtain a wavefunction 
incorporating intramonomer electronic correlation, | x(1 ),..., x(n t )) c by 
associating to each monomer not implied in a Tq-LC, a state such as |-) : 



N-2n, 

x(1),...,x(n t )) c = £ Y, aN ' 2n '' ibjD y(l)- D y(j)\ x ( 1 ),-^(n t )}. (22) 
j=0 {y(1),..,y(j)} 
y(i)*x(1),...,x(n t ) 



It is straightforward to verify that these states are eigenvectors of H PPP in the 
subspace spanned by the collective excitations with a fixed number 
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Inter monomer electronic derealization. From these \x(1 ),...,x(n t )) c 

we build new collective excitations including intramonomer electronic 
correlation, noted \n t ) c , f or n t 7cy~-LCs : 



K n 'Y 1/2 J>w,...,*r*;> c = I'^- 2M '-v A /cp^|{ J -,nj) 

{*(J),..,*(« t )} 



N—2n t 



j=0 



(23) 



The energies associated to these collective excitations from \0) c are given by 



E c {n t ) = n t {E(1)-2e_) = n t (E(1)-2e c ). 



(24) 



These collective excitations interact through the following terms 



lC t{ n t y C { n t + 1 \ H PPp\ n t) C = 



(n t +1) 



(N-2n t )(N-2n t - 1) 
N-n f 



-\1/2 



a 2 (3 s . (25) 



Because of the topological constraint typical of the Tcf-LCs, the 
mathematical problem which is governed by the interaction term (25) and 
the energy term (24) is not analytically solvable, contrary to the problem of 
the intramonomer electronic correlation. To simplify further we perform an 
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additional approximation in which the interaction term I c t is replaced by a 
simplest one analogous to (19) 

J c t (n t ) = [(n t + 1)((N-1)/3-n t )f 2 43a 2 P s . (26) 

This interaction term J c t constitutes an excellent approximation If for the 
low values of n t (Fig. 3). In the following we will check that the collective 
excitations \n t ) c with n t greater than approximately N/4 are not relevant for 

the description of the ground state, so that the approximation holds even in 
an unfavourable case. 

With the interaction term Jf the problem of intermonomer 
derealization becomes formally similar to the one of intramonomer 
electronic correlations. For the later, we have seen that the system behaves 
as N coupled independent two levels systems, the coupling being (U-V)/2. 
Here, the problem is reduced to (N-l)/3 effective independent two levels 
systems \a) and \b). These states have the energies E a = and 

3a j3 s . For the 

polymer ground state we have only to consider the lowest eigenstate in 
energy, noted . This state can be easily expressed : 

= a t \a) + b t \b) (27) 

which is associated with the energy 
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e t = ^§^- 1 -iE 1 -2e c f + 12a 4 [S? (28) 



where 



a t = i P VW " Z o and b t = i 1 ~ a t ■ (29) 



For convenience let's define N t = E(( N -1 ) / 3) where E takes the 
integer part. With the help of (27) we finally obtain an approximate 
description of the ground state : 



\GS)=Y^af'- n ^^;\n t ) c (30) 

n,=0 



associated with the energy 



E GS = Ne c + ^e t . (31) 



In order to measure the quality of our analytical approximation, we 
have performed a numerical calculation with the "exact" interaction term 
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(25). In the figure 4, we have reproduced the square of the coefficients of the 
ground state wave function on the basis of the collective excitations \n t ) c 

calculated for z = and in neglecting the coulomb terms for a long polyene. 
This parameter value corresponds to the most unfavourable choice to test 
our method. We can see that for the ground state wave function the 
gaussian approximation is indeed very good. Moreover, only the collective 
excitations \n t ) c with a low n t appear to be important. The corresponding 

energy too agrees very well (to 1%). 

Comparison with the exact result of the SSH Hamiltonian. The 

Pariser-Parr-Pople Hamiltonian admits analytical solutions if we neglect 

the electron-electron interaction term, to obtain the so called Su-Schrieffer- 

Heeger Hamiltonian [5]. In the following we compare our approximate 

result with the exact analytical which is obtained in the framework of the 

SSH Hamiltonian, versus the bond alternation parameter. In figure 5 we 

E ( ) 

have represented the ratio A{z) = — in which E SSH (z) is the exact 

e ssh{z) 

analytical result of the ground state energy. 

(i) For z = 1, our result becomes exact. This is not surprising since in 
this limiting case the monomers are totally independent. 

(ii) For z = 0, we obtain almost 92% of the exact result; this good 
agreement constitutes a surprise since our model seems rather unsuitable to 
this limit. Indeed this interesting result shows that the most important 
charge fluctuations are on the range of the intramonomer and of the nearest 
neighbours, even in the unrealistic case of a conjugated polymer in which 
all bonds would have the same length. 
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It is possible to refine further this crude description by considering 
additional GLCs; the figure 5 shows also the improvement when all the 
GLCs implying two neighbour sites are taken into account. 

For the Polyacetylene, in which the generally admitted parameter 
values give z~0,15 [1, 5], the agreement reaches about 97%. On an other 
hand, it is possible to describe other compounds as effective linear chains 
with a more pronounced bond length alternation [3]. For example the 
polyparaphenylene can be crudely associated to z~ 0,3. In such a case the 
agreement is excellent, we obtain about 99%. By looking at these good 
results, we guess that our approximation, so crude it could seem with only 
three GLCs (F-LC, D-LC and Tq-LC), keeps the essential features of the 
ground state. 

An evident flaw of our approximation at this level is the neglect of 
long-range intermonomer charge fluctuation. "Long" means here over the 
next neighbour monomer. In the SSH model, all the charge transfer LCs (at 
any distance) are degenerate. Note that, on the contrary, when electron-hole 
coulomb interaction is taken into account, the energies of the short range 
LCs are decreased and our approximation will become better. 

Finite size effect. The approximate intermonomer derealization 
energy (IDE) for the polymer fails for the smallest polyenes. For instance for 
N = 2 we obtain with this expression only 80% of the exact result. However, 
this expression can be easily improved. Indeed, for N = 2 the IDE is 
obviously given by 




(32) 
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and not by — . We can then rewrite the IDE for a polyene with N double 
3 

bonds including this correction: 



E de loc(N) = £ 2 + 




(33) 



With this new expression the IDE of any oligomer is well reproduced with 
an error always lower than 2%. This improvement of the derealization 
energy will become particularly important in the following. 

V- First excited states of (+) electron-hole symmetry - Exciton states 

In this paper we only will consider now the lowest excited states of (+) 
electron hole symmetry. They are the relevant states in order to study the 
threshold of the linear absorption. We have generated above a subspace 



configurations. For the excited states, we select an other model subspace 
based on electronic configurations which differ of | {n D ,n t }) only by localized 

perturbative areas. So, in our formalism, the excited states are composed of 
two different parts: 

(i) A local zone - called the core of the excitation - which is a local 
perturbation of the ground state system. The description of this local 



relevant for the ground state 




electronic 
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excitation requires new GLCs of (+) electron-hole symmetry. We will 
examine them below. 

(ii) Outside the core of the excitation, the dynamics of the n electrons 
remains described by F, D and Tq-LC as in the ground state. 

Let us first introduce the new GLCs on which is based the description of the 
cores of the excitations. For the states of interest here, these GLCs are 
monoexcitations of (+) electron-hole symmetry which are expressed by the 
following creation operator : 

Tct(r) = —(A\ t fl t +A+ ,5 , + A\B t +A + ,5 ,) (34) 

n\ ) 2\ n+r\ n \ n+r-l ni n\ n+r\ n-l n+rl ] v ' 

This operator creates a local mono excitation in which a hole and an 
electron are r monomers apart. These GLCs will be noted Tc^-LC. Their 
energies depend on r : 

E(0) = 2/3 d + 1 ^ (35) 

E(r>0) = 2(3 d + V + A(r). (36) 

The model subspace for excited states is spanned by the complete set of 
electronic configurations with one charge transfer excitation of (+) 
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symmetry located on sites t(n) and t(n + r), n t Tc~ -LCs located on sites 
labelled {x(i)} (i = 1,..,n t ) and n D D-LCs located on sites labelled {y(j)} 
(j = 1 ,..,n D ). Of course the same constraints as seen for the ground state have 
to be considered i.e. no spatial overlap of the different LCs can occur. A 
general expression is then 

\t(n),t(n + r);x(1),...,x(n t );y(1),...,y(n D )) = Tc+(r)Tc+ (1) ...Tc+ (rh) Dy (1) ...D+ (nD) 

(37) 

The diagonalization into this model subspace follows the previous 
approximate treatment performed for the ground state. 

The excitations (37) interact via H ppp . As for the ground state we begin 

by introducing intramonomer electronic correlations : we associate a state 
|-) (21) to every F-LC of \t(n),t(n + r); x(1 ),..., x( n t )) , to obtain 

\t( n),t(n + r); x(1 ),...,x( n t )) c = 
N-2n,-2 

£ ^a N ~ 2n >~ 2 ~ j b j D; (1) ...D; (j) \t(n),t(n + r);x(1),...,x(n t )) 

j &(l),-,y(j)} 

y( i)*t( n),t( n+r),x(1 ),..., x( n, ) 

(38) 

It is important to notice that, because of the presence of the charge transfer 
excitation, the number of F-LCs is in this case only N-2n t -2 (or 
N - 2n t -1 if the excitation is concentrated onto only one monomer, r=0). 
Indeed it is not possible to introduce D-LCs in the monomer sites occupied 
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by this perturbation. Consequently, this results in a loss of ground state 
intramonomer correlation energy which then increases the excitation 
energy [9]. 

The collective excitations (38) interact via the transfer integral in a 
similar way to the collective excitations (23) relevant for the ground state. 
However the problem is here more complex. Indeed the presence of a 7c*- 
LC divides the chain in three different sections: 

- the part on the left of the 7c* -LC which contains N L monomers, 

- the part on the right of the 7c + -LC which contains N R monomers, 

-the part inside the 7c + -LC which contains Nj monomers located 
between t(n) and t(n + r). 

We have seen before that the expression (33) is valid (with a very good 
agreement) for any chain size. The three parts of the chain can be 
independently solved in the same way than the ground state. So we obtain 
for the three parts 




(39) 



Edeloc(N|) = e 2 + 



(N,-2) 



e t , if N| > 1 



3 



(40) 



= 



else. 
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The total intermonomer derealization energy is the sum of these three 
different parts; this is a function of r, the distance between the hole and the 
electron of the excitation 7c + -LC. 

E Moc( r) = (2 + &{ r-3))e 2 + »-6-e(r-2)-e { r-3) ei m 

where @{x) is the Heavyside's function which is equal to zero for x below 
zero and one for x above zero. In every case there is a loss of IDE in the 
excited states by respect to the ground state. 

The corresponding wave function is given by 

N r 

\t(n);t(n + r)) dc = ^a Nr ' n 'b n ' \t(n);t(n + r);{n t }f (42) 
n,=0 

where N r = E(( N L + N R + N t - 3) / 3) and \t(n);t(n + r);{n t \f are collective 

excitations analogue to (23) but with local excitation on n and n + r. In this 
expression, every of the three different parts introduced by the "defect" 7c* - 
LC is independently described by a state as (30). However this approximation 
found for the ground state is relevant for the limiting case ./V— >°°. 
Consequently, the parts outside the 7c* -LC are well described in (42). On the 
contrary, the description of the part inside the 7c* -LC could become 
inappropriate when the charge transfer extends only over few monomer 
units; we have neglected these finite size effects. Moreover, in practice we 
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have neglected too the renormalization of the interaction terms due to the 
effects of intermonomer derealization. Taking into account these effects is 
straightforward without significant change in the results. 

Dressing of the excitations. Yu and co-workers dressed the particles by 
a polarization cloud following the projection technique [3, 9]. In this work, 
we dress each 7c* -LC using a perturbative treatment. In order to perform 
this, for a given 7c + -LC we consider the local configurations |7) called 
corrective local configurations (CLC), with excitation energy E lr which 
interact directly with it through tj. The effect of |7) on the energy of the 
excitation is then accounted for by a simple second order perturbation 
expansion which depends on r 

-TT ( 43 ) 

, E(r)-Ej 

The relevant CLCs give two distinct contributions to the polarization cloud: 

(i) the first kind of contribution is the dipolar one which is due to the 
long range part of the Coulomb term; 

(ii) the second kind of contribution introduces additional kinetic terms 
through the transfer integral. 

The energy of a local excited configuration, including intramonomer 
electronic correlations and intermonomer electronic derealizations, is then 
given by the following sum 
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E Dc (r) = E{r) + (N-1-0(r- 2))e c + E deloc (r) + S p (r) . 



(44) 



This energy can be expressed by respect to the ground state energy 



e(r) = E Dc (r)-E GS = E(r)-(1 + e(r-l))e c + (1 + e(r-3))e 2 

-(3+0{r-1) + 0(r-2))e t /3 + e p {r) 



(45) 



In this expression three different competing terms contribute to decrease or 
increase the excitation energy compare to the crude energy E(r). There is a 
loss of ICE and IDE and a gain of polarization energy. 

Renormalization of the interaction terms. The effective interaction 
terms between 7c* -LCs are also easily introduced by the quasi-degenerated 
second order perturbation theory. Indeed, let us note |7) and \2) two 
7c*-LCs of energies E' and E" respectively. These two Tc^-LCs interact with 
some identical CLC \l) through interaction terms noted t u and t 2 i- In these 
conditions an effective interaction exists between \1) and \2), which is 
expressed by the following simple expression 




(46) 
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This correction modifies the crude interaction term between the two LCs \1) 
and \2). 

Diagonalization. Once the energy terms (45) and the interaction terms 
have been determined, we obtain the analogue of the Rice and Gartstein's 
model for PA but with parameters directly expressible from H PPP . We have 
seen above the simple expression of V RG . We can now express the 
"correlation energy gap" with the help of the H PPP parameters: 



The excitation spectrum is then easy to calculate. Because the system is 
translationally invariant we build first, from the local excitations (42), the 
collective excitations characterised by the wavenumbers k: 



The second step consists in diagonalizing the following tridiagonal 
matrix 



U RG = lim(e(oo)-e(0)) = V+e p (oo)- £p (0) = V- 



2T(1) 2 
(U-V) 



(47) 




(48) 



n 
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(k,r\H PPP \k',r') = S w {8 rr ,e{r) + 8 r >, r±1 I r , r±1 {k)} (49) 

where 5 is the Kronecker's index and I rtr +i{k) the interaction term between 
\k,r) and \k,r±1). 

Because the interactions are local, only the collective excitations with 
the same wave number k may interact. Here we have to consider only the 
k=0 subspace in which the whole oscillator strength is concentrated. 

Size-consistency of our procedure 

Excitonic calculations are usually performed at the SCF level via a 
configuration interaction (CI) in the subspace of monoexcitations [6]. It is 
well known that size-consistency of such a procedure is not insured and 
these methods fail for instance when biexcited configurations are 
considered [7]. This is the case for the study of 2Ag states of some conjugated 
polymers [1,2]. The reason is indeed that it is impossible to treat the ground 
state and the excited states on an equal footing: doubly excited configurations 
introduce some electronic correlation in the ground state; on the contrary, 
the excited states (1B U , 2Ag) stay at an uncorrelated level. Consequently, the 
excitation energies diverge when the system size increases. This is the case 
for any incomplete CI procedure not restricted to monoexcited 
determinants. 

In our procedure, the electronic system is described at a local scale. 
The total energy is an extensive quantity and the only differences between 
the ground state and the excited configurations (equations 30 and 42) are 
localized on the core of the excitation. The other part of (42) is treated exactly 
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in the same way as for the ground state. Consequently, our procedure, as the 
one of the ref [3], is size-consistent. The introduction of higher excitation 
processes becomes straigthforward: it would suffice to consider some other 
GLCs for building the corresponding core excitations. 

Comparison with the Yu et al. results. In view of describing the first 
excited states of the PPV, Yu et al. have considered an effective Pariser-Parr- 
Pople Hamiltonian [3]. Assuming that the benzene rings just affect the 
electron transfer between double bonds, they have considered the PPV as a 
PA with a slightly more pronounced dimerization; they took z = 0.19. They 
performed a calculation of the excitation spectrum in two different steps : 

- first, they determined the band structure, taking into account the 
effect of the electronic correlations by applying the projection technique of 
Becker and Fulde [8, 9] ; 

- second, using the correlated band structure they performed an 
excitonic calculation following Abe et al. [6]. 

With z = 0.19, U = 3t , V = t and t = 2eV, the lowest singlet exciton 
of (+) electron-hole symmetry appears at E exc = 2.42 eV and the edge of the 
conduction band appears at E gap =3.58eV. These calculated values agree 

well with the experimental ones. With our method the calculation of the 
excitation spectrum for the same parameters gives E exc = 3.30eV and 
Egap = 3.55 eV. The agreement for the band gap is excellent, contrary to the 
exciton energy value, and this pinpoints a fundamental difference between 
the two models. We have to note before going further that the value of E exc 
is calculated on a finite polyene with N=75 double bonds. However the 
convergence of the calculation with the system size is fast and the finite size 
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does not affect noticeably the final results. The discrepancy between the two 
calculations can rather be attributed to the renormalization procedures 
which are radically different. In the method of Yu et al. the hole and the 
electron are renormalized by the projection technique independently each 
other [9]. Therefore the hole and electron polarization clouds are invariant 
by respect to the separation distance r between them. The energies e(r) of 
the charge transfer excitation vary only by the effect of the attractive term 
between the hole and the electron (s(r) < e(r') if r<r). This 
approximation is justified whenever the description in term of bands is 
sufficient; then electron and hole are far away on the average. In the 
excitonic states, they are constraint to stay close to each other and, moreover, 
it is known that the exciton radius is rather small in conjugated polymers. 
In this extreme case the polarization clouds of the two quasi-particles may 
interact notably and the approximation of Yu et al. becomes questionable. In 
our method we adopt a completely different procedure. Indeed we have 
renormalized the local excitation 7c* -LC in which a hole and an electron 
are separated by r monomer units. Then the hole and the electron are 
renormalized together and the correction £ p (r) depends on the distance r 

between them. When the excitation is a band to band excitation, the two 
procedures are equivalent because we can consider that the two quasi- 
particles are (principally) far from each other. So we obtain approximately 
the same value for E gap . On the contrary, for an excitonic state, the two 

procedures are different and so give different values for E exc . Yu et al. 
should have overestimated the gain of polarization energy. Indeed, when 
the two particles are very close, we can roughly guess that in the 
intermediate zone the polarization gain is counted twice. 
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In view to confirm these conclusions, it is easy to adapt our formalism 
to the hypothesis of Yu et ah. In this case we write the excitation energy of a 
charge transfer excitation of radius r £ (r ) as 



e y (r>0) = e(oo) + A(r) 



(50) 



e y (0) = e(oo) + 



U-3V 
2 



(51) 



where £(°°) is calculated following the expression (45). We have yet to 

neglect any difference between the interaction terms in the matrix (49) and 
to take I r r±1 (k) = ^(k), the interaction term between charge transfers of 

infinite radius. With these new values, we obtain E g = 3.55eV and 

E u = 2.47 eV . This excellent agreement completely confirms the origin we 

have assigned to the discrepancy. 

Exciton versus conduction band absorption. Let us finally introduce the 
interaction of an electric field E(r,t) with the 71 electrons via the dipolar 
approximation 



H int =erE{r,t) 



(52) 
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where e is the electron charge. This approximation is justified as long as the 
wavelength X of the electric field is greater than the characteristic length a; 
this is effectively the case in the visible range. 

The states of different electron-hole symmetry are coupled through 
H jnt . The linear absorption spectrum at low energy is entirely determined by 
the monoexcited states that we have presented in this section. In figure 6 we 
have represented the component along the molecular axis of the calculated 
oscillator strength for z = 0.19, U = 3t , V = t and t =2eV. As it is 
observed experimentally for the PDA, the major part of the oscillator 
strength is concentrated in the excitonic peak. In counterpart, the 
continuum of electron-hole becomes quasi-invisible. Results with other 
parameters show that, larger the binding energy of the exciton is, smaller the 
range of the prominent charge transfer states and the higher the intensity of 
the associated transition will be. Similar results have been found very 
recently but in a strong correlation approximation not really appropriate for 
conjugated polymers [22]. 

VI- Conclusion 

Intramonomer transfer integrals are larger than intermonomer ones 
in conjugated polymers. In this work, we have taken intentionally 
advantage of this characteristic feature of the conjugated polymers to build 
the electronic configurations of the polymer from the monomer orbitals. 
We have then diagonalized the Pariser-Parr-Pople Hamiltonian on a 
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reduced Hilbert space. By doing this, from PPP Hamiltonian, we obtain for 
the excited states calculations a simple effective molecular exciton method. 

With the adopted local description, each polymer electronic 
configuration is a distinct combination of several local electronic 
configurations. For the ground state, we have decided to retain only three 
kinds of local configurations, the so-called generative local configurations : 
F, D and Tq-LC. The first GLC represents a monomer in its ground state; 
the second and the third ones introduce intramonomer correlation and 
intermonomer derealization, respectively. These three GLCs permit to 
build a configuration subspace, on which we give an approximate analytical 
solution of Hppp. This solution is in good agreement with the exact result of 
the mean field approximation of the Su-Schrieffer-Heeger Hamiltonian. 

Excited states of low energy are then dealt with as local perturbations of 
the ground state so obtained. The composition of the perturbative zone (the 
core of the excitation) depends on the nature of the considered excitation. In 
this paper we only consider the lowest excitation of (+) electron-hole 
symmetry generated by 7c* -LC, in which one electron is transferred from a 
monomer n to the monomer (n+r). Furthermore, these 7c* -LCs are dressed 
by a perturbative treatment. The resulting excited state energies are 
determined by the losses of correlation and derealization by respect to the 
ground state, as well as by electron-hole attraction and the counteracting 
polarization energy. The results obtained in this manner possess the 
characteristics expected for conjugated polymers: the first excitations are 
excitons of short radius and the oscillator strength is essentially concentrated 
into these excitonic transitions. Our calculations moreover show the 
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necessity to properly consider the interaction between the two quasi-particles 
which constitute the exciton, especially when its radius is small. 

For simplicity, this approach has been applied here to the neutral trans- 
polyacetylene. In reason of their topology, other compounds need a more 
complex description of monomers and a larger number of GLCs. The 
polymer excited states will stem from the various possible excitations of a 
monomer. The corresponding collective excitations depend then on their 
mutual coupling, particularly on the relevant charge transfer integrals, 
whose some can be vanishing small. Several distinct excitonic states are 
then possible, besides excitations remaining mainly localized on monomers. 
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Figures captions 



Figure 1 

Representation of a polyene by N two-levels systems 
Figure 2 

PPP interactions expressed on the basis of monomer self consistent orbitals 
Figure 3 

Variation of the interaction term between the collective excitation \n t ) c et 

I \ C 

|n t +1) versus the number of charge-transfer excitations n t for a N=200 
polyene. exact term (25) ; approximated term (26). 

Figure 4 

Weight of the collective excitations \n t ) c in the ground state of a chain of 168 

monomers (without alternation nor electron-electron interaction). 

Open circles: numerical result 

Solid circles: gaussian approximation 

Figure 5 

Ratio of the ground state energy (31) to the exact analytical SSH solution 

versus the bond alternation parameter z 

Solid circles: using three GLCs (F-LC, D-LC and TcJ -LC) 

Open circles: using in addition all the GLCs implying next-neighbours 

Figure 6 

Component of the linear absorption along the molecular axis, in arbitrary 
units. 

The exciton peak and the threshold of the conduction band have been 
calculated for z=0.19, tO=2eV, U=3t0, V=tO 
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